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ABSTRACT: We perform state-of-the-art calculations for a prototypical dye sensitized so¬ 
lar cell: catechol on rutile TiO2(110). Catechol is often used as an anchoring group for 
larger more complex organic and inorganic dyes on Ti02 and forms a type II heterojunctions 
on TiO2(110). In particular, we compare quasiparticle (QP) GqWo with hybrid exchange 
correlation functional (HSE) density functional theory (DFT) calculations for the catechol- 
rutile Ti02(l 10 ) interface. In so doing, we provide a theoretical interpretation of ultraviolet 
photoemission spectroscopy (UPS) and inverse photoemission spectroscopy (IPES) experi¬ 
ments for this prototypical system. Specifically, we demonstrate that the position, presence, 
and intensity of peaks associated with catechol’s HOMO, intermolecular OH -0 bonds, and 
interfacial hydrogen bonds to the surface bridging O atoms ( 0 /,^H-C and O^^H- 0 ) may 
be used to fingerprint deprotonation of catechoTs OH anchoring groups. Furthermore, our 
results suggest deprotonation of these groups, while being nearly isoenergetic at high cov¬ 
erages, may significantly increase the photovoltaic efficiency of catechol-TiO2(110) inter¬ 
faces. 
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1. INTRODUCTION 

Catechol on TiO2(110) is a prototypical system for modelling in¬ 
dustrially relevant dye sensitized solar cells. This is for two rea¬ 
sons. ( 1 ) Catechol is often used as an anchoring group ^ for larger 
more complex organic and inorganic dyes. ( 2 ) Catechol forms type 
II heterojunctions, ^ where the dye’s highest occupied molecular or¬ 
bital (HOMO) is a gap state, while the lowest unoccupied molecu¬ 
lar orbital (LUMO) is above the conduction band minimum (CBM) 
of the substrate, e.g., Ti02. The level alignment for this system 
has previously been studied experimentally via ultraviolet photoe¬ 
mission spectroscopy (UPS)^’^ and inverse photoemission spec¬ 
troscopy (IPES),^ and thus provides an excellent benchmark type 
II system to be studied theoretically using many-body quasiparticle 
(QP) techniques, such as GoWq. For such type II systems, the 
HOMO-LUMO and HOMO-CBM separations determine the onset 
of the absorption spectra. In such photovoltaic devices, the dye’s 
role is to reduce the onset of the absorption spectra to sub-band gap 
energies, and maximize the overlap with the solar spectrum. 

The exact structure of the catechol-rutile TiO2(110) interface is 
difficult to identify and control experimentally. This is because 
for hydroxylated molecules such as catechol, a complex network 
of interfacial and intermolecular bonds is formed upon adsorption. 
Moreover, catechoTs anchoring OH groups can be either fully dis¬ 
sociated, partially dissociated, or intact, with each catechol over¬ 
layer having similar adsorption energies. Still, dissociation plays a 
significant role in determining the level alignment and, hence, 
the photovoltaic efficiency of the interface. 

Generally, the HOMO moves to higher energy upon deprotona¬ 
tion for hydroxylated molecules, e.g., CH3OH on TiO2(110) and 


H2O on rutile TiO2(110) and anatase TiO2(101). For this rea¬ 
son, the homo’s energy provides a fingerprint of the interface’s 
structure. This is augmented by the presence of distinct identi¬ 
fiable levels associated with intermolecular and interfacial OH -0 
hydrogen bonds. Specifically, the hydrogenated bridging O atoms 
(HO^^) resulting from interfacial deprotonation of the anchoring 
groups are consistently at ~ 10 eV below the Fermi level e/7 on 

To describe both highly hybridizedand localized molec- 
ular interfacial levels, one requires a correct description of the 
anisotropic screening, i.e., electron-electron interaction, at the in¬ 
terface. This is clearly seen from the level alignment of H2O on 
TiO2(110).‘® 

On the one hand, H20’s Ibi level becomes highly hybridized 
with the substrate upon adsorption on TiO2(110). The interfacial 
level alignment of these highly hybridized levels with the substrate 
VBM is already described by density functional theory (DFT). This 
is because the screening of these levels is essentially the dielectric 
constant of the substrate. At the GqWo level, the DFT Kohn-Sham 
(KS) eigenenergies are shifted by A, the difference between the QP 
self-energy (E = iGoWo) and the xc potential (Vxc), up to a normal¬ 
ization factor Z, i.e., A = Z(E - Vxc)- For example, A 0.1 and 
0.7 eV for all levels of the TiO2(110) substrate from GqWo based 
on DFT calculations employing a generalized gradient approxima- 

tion (PBE) and a hybrid range separated (HSE) exchange and 
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correlation (xc)-functional, respectively. 

On the other hand, H20’s 3 ai and lb2 levels remain mostly 
localized on the molecule upon adsorption on TiO2(110).^^ The 
alignment of these localized levels, whose screening is signifi¬ 
cantly different from that of the substrate, is poorly described at 
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the DFT level. Not even HSE provides an accurate description 
of the interfacial level alignment for localized levels, e.g., H 20 ’s 
3ai and lb 2 and CH 30 H’s a" on TiO 2 ( 110 ). This is because 
HSE DET calculations effectively perform a constant static screen¬ 
ing of the exchange term, i.e., the fraction of Hartree-Fock exact- 
exchange included a = 0.25 is effectively an inverse dielectric con¬ 
stant a For this reason, HSE performs well for systems 

with a homogeneous screening and Soo ~ 4. 

To correctly account for differences in screening between the 
molecular layer and substrate, one may use many-body QP tech¬ 
niques, such as Go Wo 2^ In such methods, the spatial dependence 
of the screening is included explicitly. In fact, for occupied lev¬ 
els, A is linearly dependent on the fraction of the wave function’s 
density within the molecular layer This means by just 

rigidly shifting all the levels, one cannot describe the alignment of 
occupied levels with significant density outside the substrate. 

For type I interfaces, i.e., H 2 O or CH 3 OH on Ti 02 (l 10), 1^46,20,23 
we have previously demonstrated that Go Wo provides an accurate 
alignment for both localized and highly hybridized levels. In each 
case. Go Wo shifts the localized levels to stronger binding, into 
quantitative agreement with UPS experiments. 

In this study, we compare the projected density of states (PDOS) 
onto catechol obtained from QP PBE Go Wo and HSE DFT calcu¬ 
lations with the measured UPS and IPES spectra for the catechol- 
rutile TiO 2 ( 110 ) interface. In so doing, this study provides a com¬ 
plete state-of-the-art computational description of the simplest ex¬ 
perimentally relevant type II interface. Based on our analysis of 
the PDOS, we are able to suggest the most likely structure of the 
catechol overlayer measured in UPS experiments. In fact, our re¬ 
sults suggest the degree of catechol dissociation may differ with the 
experimental conditions employed. Our results suggest fully de- 
protonating the anchor groups of the overlayer should lead to an 
increased efficiency of the photovoltaic device. 

2. METHODOLOGY 

Our Go Wo calculations have been performed using vasp within 

the projector augmented wave (PAW) scheme. The Go Wo calcu¬ 
lations are based on KS wave functions and eigenenergies from 

DFT obtained using PBE.^^ DFT calculations employing the 
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HSEO 6 variant of the HSE xc-functional have been carried out 
for comparison with PBE Go Wo calculations. 

In the Go Wo approach, the contribution to the KS eigenval¬ 
ues from the xc-potential Vxc is replaced by the QP self energy 
E = iGoWo in a single step, where Gq is the Green’s function and 
Wo is the screening based on the KS wave functions and eigen¬ 
values. The dielectric function is obtained from linear response 
time-dependent (TD) DFT within the random phase approximation 
(RPA), including local field effects. From Go Wo one obtains first- 
order QP corrections A to the KS eigenvalues, but retains the KS 
wave functions. 
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The geometries have been fully relaxed using the PBE xc- 
functional, with all forces < 0.02 eV/A. We employ a plane-wave 
energy cutoff of 445 eV, an electronic temperature kBT ^ 0.2 eV 
with all energies extrapolated to T ^ 0 K, and a PAW pseudopo¬ 
tential for Ti which includes the 3and 3/?^ semi-core levels. The 
calculations have been performed spin unpolarized. All unit cells 
contain a four layer TiO 2 ( 110 ) slab, employ the measured lattice 
parameters of bulk rutile Ti 02 {a = 4.5941 A, c = 2.958 A),^^ and 
include at least 27 A of vacuum between repeated images. In each 
case, equivalent catechol overlayers are adsorbed on both sides 
of the slab. We employ F centered k-point meshes with densities 
Ak < 0.25 A“2 approximately 9Va unoccupied bands per atom, i.e. 
including all levels up to 30 eV above the VBM, an energy cutoff 


of 80 eV for the number of G-vectors, and a sampling of 80 fre¬ 
quency points for the dielectric function. The GqWo parameters are 
consistent with those previously used for describing both rutile and 
anatase Ti 02 bulk, rutile TiO 2 ( 110 ) and anatase TiO 2 ( 101 ) clean 
surfaces, and their interfaces. These parameters have been 

shown to provide accurate descriptions of bulk rutile and anatase 
optical absorption spectra, and both clean surface and interfacial 
level alignment. 

The adsorption energy Eads of catechol on Ti coordinately unsat¬ 
urated (TicM^) sites of a TiO 2 ( 110 ) surface is given by 

E .„.. g[»C.«chol.T.O,(110)]-E[-nO,(110)]_^ 

n 

(1) 

where n is the number of adsorbed catechol molecules in the 
supercell, and E[fzCatechol-r TiO 2 ( 110 )], E[TiO 2 ( 110 )], and 

[Catechol] are the total energies of the covered and clean sur¬ 
faces and gas phase catechol molecule, respectively. For catechol 
in the gas phase, we find the most stable conformation has an 
intramolecular hydrogen bond, i.e., Cg symmetry. 

Scanning tunneling microscopy (STM) simulations have been 
performed using the Tersoff-Hamann approximation. In this ap¬ 
proach, the current / at a position r is given by 

r Ep + U 

/(r) = C / p{Y,s)ds^CUp{Y,SF)^CUp{Y,SF + U), (2) 

Jef 

where C is a prefactor which depends on the DOS, surface work 
function, and tip radius. U is the potential of the sample relative to 
the tip in the experiment, i.e., the applied potential relative to the 
experimental Fermi level sf- p(r,^) is the local DOS, given by 


P(r,e) = 2‘^(®nk“®) 

nk 


Nk 


«y;exp 

nk 


(^nk |^«k(*")| 

V keT ) \ 


( 3 ) 

where is the GoWo eigenvalue and is the KS wave function 
of level n at k-point k, 0.2 eV is the electronic temperature of 
the calculation, and N\^ is the weight of k-point k. To emphasize any 
dependence of /(r) on the applied bias U, we have used throughout 
/(r) CUp(y,sf + U). It should be noted, however, that similar 
results are obtained by integrating p over the bias window [sf,sf + 
C/] (see Figure SI in Supporting Information). Herein, p is plotted 
at an isosurface value of5xl0 ej A . This is somewhat greater 
than its maximum far from the surface, i.e.. 


min(maxp(v,y,z;e/ 7 -r C/) I < 5 X 10 ^ ^/A^. (4) 

^ V / 


This ensures p is defined at this isosurface value throughout the 
surface plane. 

Experimental spectra are typically referred to the Eermi level, 
sf, which is pinned ~ 0.1 eV below the CBM for mildly re- 

duced Ti 02 . Using the electronic band gap for rutile Ti 02 of 
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3.3 ±0.5 eV obtained from electron spectroscopy measurements, 
the experimental VBM energy relative to the Eermi level is £vbm ^ 
0.1 -3.3 -3.2 eV.^^ Since the VBM is the most reliable theo- 
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retical energy reference, we subtract ^vbm ^ “3.2 eV from the 
measured UPS and IPES spectra to align with the calculated GqWo 
DOS and PDOS, and vice versa. 

We align the DOS and PDOS with respect to the deepest Ti semi¬ 
core level (3s^). This allows a direct comparison between spectra 
for half {ViT) ViT)), mixed (D ViD), and fully (D D) dissociated cat¬ 
echol overlayers. In each case, the highest occupied levels belong 
to the catechol overlayer. As we are interested in seeing the dif¬ 
ference between the HOMO position and the TiO 2 ( 110 ) VBM for 
each catechol overlayer, the highest occupied level is not a good 
reference. 

Moreover, the catechol anchor groups form Tic^^-O bonds, mak- 
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ing it impractical to separate the O Ipn surface levels from those 
of the molecule. This makes it difficult to identify a highest oc¬ 
cupied level with purely surface contributions, which can be asso¬ 
ciated with the clean surface’s VBM. For this reason, one should 
use the Ti levels as reference, e.g., occupied semi-core Ti 3s^ or 
unoccupied Ti 3d CBM levels. Here, we align with respect to the 
deepest semi-core Ti 3s levels. In this way we obtain a consistent 
alignment relative to the surface levels for each configurations. We 
find this is effectively equivalent to aligning relative to the Ti 02 
CBM, i.e., Ti 3d levels. 

Using the semi-core levels, we remove differences in work func¬ 
tion between surfaces, which would be present if the vacuum level 
were used as a reference. Here, we take the energy of the VBM 
for the clean surface ^vbm relative to its deepest Ti 3s^ semi-core 
level as our final reference relative to the catechol interface’s 
deepest Ti 3s^ level More precisely, ^vbm = ^vbm + 

^TisB, Figures 3 and 4. 

However, for reduced TiO 2 -jc( 110 ), where Ti 3d levels are oc¬ 
cupied, all the Ti levels are consequently upshifted compared to 
stoichiometric TiO 2 ( 110 ). This makes the Ti levels a poor refer¬ 
ence for comparison between such systems. For this reason, we use 
the VBM as an energy reference for 1 ML H@0/;^ on TiO 2 ( 110 ), 
as it is a type I interface. 

3. RESULTS AND DISCUSSION 

Catechol consists of a benzene ring with two adjacent anchoring 
OH groups. It has been previously shown, both theoretically and 
experimentally, that catechol adsorbs on the Tic^^ sites of the rutile 
Ti 02 (l 10) surface via the OH anchor groups in a bidentate config¬ 
uration. 

At low coverage, catechol preferentially adsorbs upright (6 ^ 
86 °) on the surface, parallel to the [ 001 ] c-axis,^ with both anchor¬ 
ing groups deprotonated, i.e., fully dissociated, with an accompa¬ 
nying charge of -0.4^ transferred to the nearest O^r atom of the 
TiO 2 ( 110 ) surface. As the coverage increases (~ 2/3 ML), cate¬ 
chol tilts (6 ~ 67°) toward the surface, with two interfacial 0/,^H-0 
bonds (cf. Table 1). 

Table 1. Adsorption Energies Eads in Electronvolts per Molecule for V 2 
ML 1x4, 2/3 ML 1x3, and 1 ML 1x4 Intact (I), Half ( 1 / 2 D) and Fully (D) 
Dissociated Catechol Overlayers on TiO2(110) with Tilting Angle 6 in 
Degrees and Number of Interfacial (OH-O) and Intermolecular (OH- 
and O^^H-O) Bonds per Unit Cell 

structure 6 OH-O OH-0^^ O^^H-C O^^H-0 

(ML) (°) (bonds/unit cell) (eV) 


1/2 D 

86 

0 

0 

0 

2 

-0.792 

I 

87 

0 

0 

0 

0 

-0.249 

I 

61 

0 

2 

0 

0 

-0.595 

'/2D 

80 

0 

1 

0 

1 

-0.620 

% '/2D 

63 

0 

1 

0 

1 

-0.676 

D 

78 

0 

0 

2 

0 

-0.676 

D 

87 

0 

0 

0 

2 

-0.705 

D 

67 

0 

0 

0 

2 

-0.748 

'/2D ViD 

56, 49 

2 

0 

2 

0 

-0.614 

ViD ViD 

55,48 

1 

1 

0 

2 

-0.600 

ViD ViD 

54, 46 

2 

0 

0 

2 

-0.598 

^ D '/2D 

56, 48 

1 

0 

1 

2 

-0.685 

D '/2D 

56, 48 

1 

0 

2 

1 

-0.653 

DD 

56, 49 

0 

0 

2 

2 

-0.652 

At high coverage (1 ML), where catechol molecules adsorb on 
every Tic^^ site, they are forced to tilt in alternating directions due 

to steric hindrance. This 

gives 

rise to the 1 x 

4 catechol overlayers 

seen via STM ^ 

and shown in 

Figure 1. 

Further, at this 

coverage. 


besides interfacial hydrogen bonds to the surface HO/,^, the over¬ 
layer, when not fully dissociated, is stabilized by intermolecular 
hydrogen bonds between neighboring catechol molecules. For this 
reason, the half, mixed, and fully dissociated catechol adsorption 
energies are all within 0.1 eV, i.e., the accuracy of DFT (cf. Ta- 
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Figure 1. Schematics of 1 ML catechol adsorbed (a) half (/ 2 D YiD, blue), 
(b) mixed (D ViD, green), and (c) fully (D D, red) dissociated on coordi- 
nately unsaturated Ti sites (Ticus) of TiO2(110). Charge transfer of —0.4^ 
accompanying deprotonation is represented by arrows, while intermolecu¬ 
lar OH-O (gray) and interfacial O^^H-O (black) and O^^H-C (magenta) 
hydrogen bonds are denoted by dotted lines. The angle between catechol’s 
benzene ring and the surface plane, 6, is shown above. 


ble 1). These results have also been reproduced to within 25 meV 
using the real space PAW DFT code gpaw. ’ 

Overall, we find catechol binds more weakly to the surface as 
the coverage increases. This is attributable to steric hindrance, es¬ 
pecially for a 1 ML coverage. Furthermore, the binding energy at 
y 3 ML coverage is significantly stronger for deprotonated anchoring 
groups and tilted catechol molecules (cf. Table 1). This is consis¬ 
tent with previous DFT studies of catechol on TiO 2 ( 110 ).^ 

The most stable 1 ML half, mixed, and fully dissociated cate¬ 
chol overlayers listed in Table 1 are shown in Figure 1 . The 1 ML 
y 2 D ViD structure (Figure 1(a)) has two intermolecular OH-O (gray 
dashed lines) and two interfacial O^^H-C bonds (magenta dashed 
lines), the 1 ML D ViD structure (Figure 1(b)) has one intermolec¬ 
ular OH-O, one interfacial 0/,^H-C, and two interfacial 0/,^H-0 
bonds (black dashed lines), while the 1 ML D D structure (Fig¬ 
ure 1 (c)) has two interfacial O^^H-C and two interfacial O/^^H-O 
bonds. In each case, interfacial deprotonation of the anchoring OH 
groups is accompanied by a charge transfer of -0.4^ to the nearest 
Ohr atom, as depicted schematically in Figure 1 . Such intermolec¬ 
ular and interfacial hydrogen bonding combinations have also been 
reported for the 1 ML methanol-TiO 2 ( 110 ) interface. Figure 1 
shows that the two tilting directions are inequivalent, with adjacent 
catechols tilted by ^ ~ 56° and 48°. For the 0/,^H-0 interfacial hy¬ 
drogen bonds, catechol tilts away from the O/^^H moiety, while for 
the O^^H-C interfacial hydrogen bonds, catechol tilts towards the 
Oljril moiety. For the latter, the 0/,^H-C bond is mostly KOhr cr in 
character, with a minor C 2p;^ contribution from the neighboring C 
atoms. 

STM images of the 1 x 4 catechol overlayer, ^ and simulated im¬ 
ages for the y 2 D VzD, D VzD, and D D structures are shown in Fig¬ 
ure 2. In order to reproduce the line scan’s minima, i.e., effective 
height, it is necessary to perform the line scan along the [ 111 ] di¬ 
rection, as indicated by arrows in Figure 2(a-f). This is particularly 
important for the D D structure, where a line scan along the [111] 
direction would not cross the computed STM minima. Comparing 
the measured and computed line scans along the [ 111 ] direction, 
we thus conclude that the measured catechol overlayer is the mirror 
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Figure 2. Constant current (a-f) STM and (g-1) line scans at (left panels) U = 0.6 V and (right panels) U = 0.9 V from ref 1 (upper panels and brown thick 
lines) and calculated LDOS (lower panels and blue/green/red thin lines) for (a,d,g j) half (ViD ViD, blue), (b,e,h,k) mixed (D ViD, green), and (c,f,i,l) fully (D 
D, red) dissociated catechol overlayers on TiO2(110) at 5 x 10“^ e/A?. Heights in A are relative to the measured or calculated minimum. 


image of the computed structures shown in Figure 1. 

Overall, the average heights h from line scans and simulated 
STM images agree qualitatively for all three computed structures at 
both U = 0.6 Y{h^ 2.79 A, /texp ^ 2.87 A) and U = 0.9Y{h^ 2.67 
A, /texp 2.75 A). However, the difference in height between neigh¬ 
bouring molecules Ah is underestimated at C/ = 0.6 V (Ah 0.14 A, 
A/iexp ^ 0.81 A), but agrees qualitatively at C/ = 0.9 V (Ah ^ 0.32 A, 
A/iexp 0.22 A). Altogether, this indicates that the structure of the 
unoccupied levels at these biases are rather insensitive to the depro¬ 
tonation of catechol’s anchoring OH groups. This is not surprising, 
considering the structural similarity of the benzene ring orientation 
for each type of overlayer shown in Figures 1 and 2. 

In summary, STM provides direct information as to the rela¬ 
tive orientation of catechol on the surface which forms the 1x4 
overlayer on TiO 2 ( 110 ). However, STM lacks direct information 
about deprotonation of the OH anchor groups. Such information 
is important, as it determines the relative level alignment of the 
molecule’s HOMO with the substrate’s VBM. To obtain direct in¬ 
formation about the level alignment, one must compare the GqWq 
PDOS with UPS and IPES spectra. 

In Figure 3, we identify features of the GqWq PDOS which may 
be used to fingerprint the catechol overlayer’s structure. These 
features are the HOMO, HOMO-1, interfacial and inter- 

molecular HO-0 levels associated with the half (blue/red) and fully 
(green/orange) dissociated surface catechol species. Dotted lines 
connect energy levels associated with the half or fully dissociated 
catechol species. 

Although catechol’s gas phase spectrum contains many addi¬ 
tional occupied tt and cr skeletal levels, shown in gray in Fig¬ 


ure 3(a), these levels are rather insensitive to deprotonation of the 
OH anchor groups. For this reason, they are ineffective for distin¬ 
guishing between ViD V 2 D, D V 2 D, and D D catechol overlayers on 
TiO2(110). 

As a reference, the HOMO and HOMO-1 levels of gas phase 
catechol are depicted in Figure 3(a), aligned relative to the 
molecule’s HOMO. Figure 3(b), (c), and (d) show that the HOMO 
and HOMO-1 are pinned to each other, and shift up in energy 
with deprotonation of adsorbed catechol. This deprotonation of the 
OH anchor group induces a charge transfer of —0.4^ to the sub¬ 
strate. As charge is removed from the molecule, the HOMO and 
HOMO-1 are destabilized. This effect is even more pronounced at 
the GqWq level, as the molecule’s ability to screen the HOMO and 
HOMO-1 levels is also reduced as charge is transferred to the sub¬ 
strate. Consequently, the energy separation between HOMO and 
VBM may be used to distinguish between half, mixed, and fully 
dissociated catechol overlayers on TiO 2 ( 110 ). 

Another fingerprint of dissociated catechol is the presence of 
KOijr surface levels at ~ 10 eV below the experimental Fermi level. 
The HObr level is a general feature of all the interfaces formed 
from rutile TiO 2 ( 110 ) and hydroxylated molecules, e.g., H 2 O and 
CH 3 OH. As a reference, we show in Figure 3(e) both in and 
out of phase HObr cr levels for 1 ML H@0/,^, on TiO 2 ( 110 ).^^ 
This structure is equivalent to 1/2 ML dissociated H 2 O adsorbed on 
bridging O vacancies (H 20 @ 0 ^^^) of a reduced TiO 2 -i/ 4 ( 110 ) sur¬ 
face. In Figure 3(c,d) we show HO/,^ levels associated with ad¬ 
jacent 0/,^H-0 bonds, while for (b), we show O^^H-C levels. For 
this reason, the HO^^ levels in (b) have less weight on the benzene 
ring compared to (c,d). Despite the differences in reduction of the 
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Figure 3. HOMO, HOMO-1, OH-O, and HO^^ DFT orbitals and GqWq eigenvalues and DOS for catechol (a) in the gas phase (grey lines), 1 ML (b) half 
(ViD ViD, blue), (c) mixed (D YiD, green), or (d) fully dissociated (D D, red), and (e) H@0^r (grey) on TiO2(110). (b,c,d) Colors are used to differentiate 
between YiD (red/blue) and D (orange/green) catechol orbitals. Energies are relative to the VBM (^vbm) of clean TiO2(110) (left) or the experimental Fermi 
level (right). 


substrate between the systems depicted in Figures 3(b,c,d,e), the 
KOijr levels are surprisingly consistent in energy. 

As can be seen in Figure 1(a), for ViD V 2 D, the groups are 
on every other whereas in Figure l(b,c), the HO/,^-0 bonds 
for D V 2 D and D D are on adjacent For this reason, the HOtr 
levels are essentially isoenergetic in Figure 3(c,d,e), whereas in (b) 
the HObr level is more stable by ~ 0.3 eV. This is consistent with the 
observed downshift by 0.1 eV of the HObr surface levels upon re¬ 
ducing the coverage from 1 ML to ^2 ML ll@Obr on Ti02(l 10). 

A peak at ~ 12 eV below the experimental Fermi level indicates 
the presence of intermolecular OH-O hydrogen bonding within the 
overlayer. On the one hand, this may be used to fingerprint the 
presence of catechol which is not fully deprotonated. On the other 
hand, its absence suggests the catechol overlayer is fully dissoci¬ 
ated. These levels have significant cr-bonding character along the 
OH-O intermolecular hydrogen bond. This is combined with ben¬ 
zene skeleton cr orbitals. For the D ViD overlayer, the OH-O orbital 
has very little weight on the fully dissociated catechol molecule, as 
seen in Figure 3(c). In other words, the OH-O levels are mostly 
associated with V 2 D catechol molecules. 

In Figure 4, we compare the GqWq DOS and PDOS for the clean, 
ViD ViD, D ViD, and D D catechol overlayers with IPES ^ and UPS ^ 
for the (a) clean, (b) catechol covered, and their (c) difference spec¬ 
tra. In each case, there are several peaks outside the clean surface’s 
VB region, shown in gray in Figure 4. These are the peaks which 


are most easily distinguishable from the surface levels. For over¬ 
layers which are highly hybridized with the surface, e.g., H 2 O and 
H covered TiO2(110), it is difficult to disentangle surface and 
molecular levels using the difference between experimental spectra 
for the covered and clean surfaces, i.e., difference spectra. Overall, 
our Go Wo DOS and PDOS agree qualitatively with the UPS/IPES 
and AUPS/AIPES spectra for catechol on TiO2(110) from refs 1 
and 9. 

Although there is a nice alignment between the computed 
and measured spectra for unoccupied levels (Eigure 4(b)), these 
substrate levels are hybridized with the catechol overlayer (Eig¬ 
ure 4(c)). As a result, although the unoccupied Go Wo DOS, PDOS, 
IPES and AIPES levels agree, there is little to distinguish between 
the types of catechol overlayers. This is not at all surprising, as 
these levels are predominantly Ti 3d in character, and should hy¬ 
bridize equally well with ViD and D catechol. However, there is a 
noticeable increase in PDOS intensity for the D D structure. This 
is probably associated with a stronger coupling between 0-Ticus 
compared to HO-Tic^.?. This results in a raising of the Ticus atoms 
out of the Ti02(l 10) surface plane in Eigure 1(c) for the D D cate¬ 
chol overlayer. 

By comparing the UPS peaks outside the VB region to the GoWo 
DOS and PDOS, we find features suggestive of V 2 D V 2 D catechol in 
the UPS of ref 9, and D D catechol in the UPS of ref 1 . However, 
differences in detection setup and resolution between refs 9 and 1 
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Figure 4. GqWq total DOS for (a) clean TiO2(110) and (b) 1 ML catechol 
covered TiO2(110), and (c) GqWq (thick lines) and HSE DFT (thin lines) 
catechol PDOS, for half {ViD ViT), blue), mixed (D ViD, green), and fully 
(D D, red) dissociated catechol overlayers, as compared to the experimental 
UPS and IPES from ref 9 (blue circles) and UPS from ref 1 (red diamonds) 
for (a) clean Ti 02 (l 10), (b) catechol covered Ti 02 (l 10), and (c) their differ¬ 
ence spectra. Energies are relative to the VBM (^vbm) of clean Ti02(l 10). 
Filling denotes occupation. 


mean an absolute attribution of the measured spectra to V 2 D V 2 D and 
D D catechol, respectively, may be excessive. Nevertheless, com¬ 
paring the UPS catechol covered and difference spectra from ref 9 
with the GqWq spectra, we find the following three fingerprints of 
the V 2 D ViD spectra. (1) A shoulder ~ 0.5 eV above the VBM is 
suggestive of the ViD V 2 D catechol overlayer’s HOMO. (2) A peak 
at ~ -10 eV is suggestive of HOtr interfacial levels. (3) A peak at 
—12 eV is suggestive of intermolecular OH-0 hydrogen bonds. 
Performing a similar comparison to the catechol covered UPS from 
ref 1 , and the resulting difference spectra, we find the same three 
fingerprints, but of the D D spectra. (1) A well-separated peak ~ 0.8 
eV above the VBM is suggestive of the D D catechol overlayer’s 
HOMO. (2) A more intense peak at ~ -10 eV is suggestive of HO^^ 
interfacial levels. (3) A significant dip in the spectra at —12 eV 
suggests fewer intermolecular OH-0 hydrogen bonds. Note that, 
as HO^^ levels are associated with the Ti 02 (l 10 ) surface, this peak 
is absent from the PDOS, as the O^r atom is part of the surface. 

In Figure 4(c), we also compare catechol’s PBE GqWq and HSE 
DET PDOS for ViD V 2 D, D V 2 D, and D D overlayers on TiO 2 ( 110 ). 
In each case, the HSE DET PDOS for occupied levels yields ~ 1 
eV weaker binding energies than PBE GqWq. This upshift of the 
occupied molecular levels with HSE DET is consistent with our 
previous results for the localized CH 30 H’s a" and H 20 ’s 3ai and 
lb 2 levels on TiO 2 ( 110 ). On the other hand, for the unoccu¬ 
pied levels, catechol’s PDOS from HSE DET and PBE GqWq are 
consistent with each other, as the unoccupied molecular levels are 
highly hybridized with the substrate. This reinforces the finding 
that HSE fails to provide an accurate description of the interfacial 
level alignment for localized molecular levels. 


4. CONCLUSIONS 

The structure of catechol overlayers on TiO 2 ( 110 ) is characterized 
by a complex network of interfacial and intermolecular hydrogen 
bonds. It is difficult to precisely establish the detailed structure of 
the catechol overlayer based solely on STM experiments. This is 
because simply deprotonating catechol’s anchoring groups, while 
nearly isoenergetic, strongly affects the position of the HOMO in 
this type II interface. Hence, the extent of catechol’s deprotonation 
on the surface determines the interface’s photovoltaic efficiency. 

We combine GqWq level alignment with UPS measurements to 
identify the fingerprints of half, mixed, or fully dissociated cate¬ 
chol overlayers on TiO 2 ( 110 ). Only QP techniques, such as GqWq, 
are sufficiently accurate to robustly predict an ~ 0.5 eV energy dif¬ 
ference between catechol’s HOMO position in the half and fully 
dissociated catechol overlayers. Moreover, besides the HOMO po¬ 
sition, which indicates the extent of deprotonation of catechol’s OH 
anchoring groups, the absence of a peak at —12 eV is indicative 
of a lack of intermolecular OH-0 bonds. Likewise, the presence 
of a peak at —10 eV suggests the presence of KOjjr groups on 
the surface, which are formed upon deprotonation of catechol’s OH 
anchoring groups. 

This distinct peak has also been observed for H 2 O dissociated on 
bridging O vacancies, i.e., H20@0^^^, of a reduced TiO2-jc(110) 
surface. ^4,16,18,19 previously been considered a fin¬ 

gerprint of dissociated H20@0^^^. However, we have shown this 
KOhr peak at —10 eV may also arise from deprotonation of hy- 
droxylated molecules on Ticus- Eurther, the position of the HO^^ 
peak is rather insensitive to the degree of substrate reduction. In 
UPS experiments, one should find this peak moves to slightly 
weaker binding energies as the coverage of HO^^ is increased. 

While the energy of the unoccupied levels is rather insensitive to 
deprotonation of catechol’s anchoring groups, the overlap of these 
substrate levels with those of the catechol overlayer increases with 
deprotonation. Combined with the destabilization of catechol’s 
HOMO with deprotonation, this suggests fully deprotonated cat¬ 
echol overlayers should have the greatest photovoltaic efficiency. 
This work provides a road map for future studies of catechol’s op¬ 
tical absorption based on the Bethe-Salpeter equation, and cat¬ 
echol’s subsequent rate of charge transport through the Ti 02 sub- 
strate using non-equilibrium Green’s function methods. ’ 
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